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Abstract 

o . 

A brane scenario is envisaged in which the m-dimensional bulk is endowed with a Gauss- 
es [ Bonnet term and localization of matter on the brane is achieved by means of a confining potential. 

The resulting Friedmann equations on the brane are modified by various extra terms that may be 
interpreted as the X-matter, providing a possible phenomenological explanation for the accelerated 
expansion of the universe. The age of the universe in this scenario is studied and shown to be 
t-H | consistent with the present observational data. 

7— I . 

1 Introduction 

Recent years have been witnessing a phenomenal interest in the possibility that our observable four- 
dimensional (4D) universe may be viewed as a brane hypersurface embedded in a higher dimensional 
bulk space. Physical matter fields are confined to this hypersurface, while gravity can propagate in 
the higher dimensional space-time as well as on the brane. The most popular model in the context 
of brane world theory is that proposed by Randall and Sundrum (RS). In the so-called RSI model 
[1], they proposed a mechanism to solve the hierarchy problem with two branes, while in the RSII 
model [2], they considered a single brane with a positive tension, where 4D Newtonian gravity is 
recovered at low energies even if the extra dimension is not compact. This mechanism provides us 
an alternative to compactification of extra dimensions. The cosmological evolution of such a brane 
^ ■ universe has been extensively investigated and interesting effects such as dark radiation or quadratic 

density term in Friedmann equations have been found [3]. 

There is a general belief that Einstein gravity is the low-energy limit of a quantum theory of 
gravity which is still unknown. A promising candidate is the string theory which suggests that in 
order to have a ghost-free action, quadratic curvature corrections to the Einstein-Hilbert action must 
be proportional to the Gauss-Bonnet (GB) term [4]. This term also plays a fundamental role in Chern- 
Simons gravitational theories [5]. The appearance of higher derivative gravitational terms can also be 
seen in the renormalization of quantum field theories in curved space-times [6]. From a geometric point 
of view, the combination of the Einstein-Hilbert and Gauss-Bonnet term constitutes, for 5D space- 
times, the most general Lagrangian producing second-order field equations [7]. In four dimensions, 
the Gauss-Bonnet combination reduces to a total divergence and is dynamically irrelevant. These 
facts provide a strong motivation for the study of braneworld theories which include a Gauss-Bonnet 
term. Interestingly, it has been shown that the graviton zero mode is as localized at low energies in 
the Gauss-Bonnet brane system as in the RSII model [8]- [10] and that the correction of the Newton's 
law becomes less pronounced by including the Gauss-Bonnet term [11]. 

As the discussion above suggests, it comes as no surprise that brane theories with a Gauss- 
Bonnet term have been generating so much interest over the recent past [12]. One of the results 
of these investigations is that the Gauss-Bonnet correction to the Friedmann equations has been 
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dominant in the early universe, having the form H 2 oc p 2 ^ 3 at high energies. This term arises from 
the imposition of the generalized Israel junction conditions [13]. However, it has been argued that 
such junction conditions may not be unique. Indeed, other forms of junction conditions exist, so that 
different conditions may lead to different physical results [14]. Furthermore, these conditions cannot 
be used when more than one non-compact extra dimension is involved. To avoid such concerns, 
an interesting higher-dimensional model was introduced in [15] where particles are trapped on a 4- 
dimensional hypersurface by the action of the confining potential V. In [16], the dynamics of test 
particles confined to a brane by the action of such potential at the classical and quantum levels were 
studied and the effects of small perturbations along the extra dimensions investigated. Within the 
classical limits, test particles remain stable under small perturbations and the effects of the extra 
dimensions are not felt by them, making them undetectable in this way. The quantum fluctuations 
of the brane cause the mass of a test particle to become quantized and, interestingly, the Yang-Mills 
fields appear as quantum effects. Also, in [17], a braneworld model was studied in which matter 
is confined to the brane through the action of such a potential, rendering the use of any junction 
condition unnecessary and predicting a geometrical explanation for the accelerated expansion of the 
universe. 

In brane theories the covariant Einstein equations are derived by projecting the bulk equations 
onto the brane. This was first done in [18] where the Gauss-Codazzi equations together with Israel 
junction conditions were used to obtain the Einstein field equations on the 3-brane. The same 
procedure was subsequently used in [19] to obtain the field equations on the brane in the presence of 
the Gauss-Bonnet term where to confine matter on the brane, the authors used a delta-function in 
the energy-momentum tensor. A feature of these works is that the effective field equations are too 
complicated to be physically transparent because the extrinsic curvature cannot be expressed solely 
by the energy-momentum tensor on the brane. Still, they are important in assessing the effects of 
the higher order terms appearing in the action. It would therefore be interesting to investigate the 
effects of such higher order terms when a confining potential, rather than junction conditions, is used 
to localize matter on the brane. 

In this paper, we consider an m-dimensional bulk space in the presence of the Gauss-Bonnet 
term without imposing the Z2 symmetry [20]. As mentioned above, to localize matter on the brane, a 
confining potential is used rather than a delta-function in the energy-momentum tensor. The resulting 
Friedmann equation is modified by the appearance of extra terms that behave like the X-matter; the 
phenomenological model proposed to fit the data explaining accelerated expansion of the universe. 
Within the framework of our model, we also calculate the age of the universe and show that it is 
consistent with the bounds set by the modern observational data. We should emphasize here that 
there is a difference between the model presented in this work and models introduced in [9, 10] in 
that in the latter no mechanism is introduced to account for the confinement of matter to the brane. 

2 Geometrical considerations 

Consider the background manifold V m -\ isometrically embedded in a pseudo-Riemannian manifold 
V m by the map y : V m -\ — ► V rn such that 

G AB yp B v = g^ g AB y A ^ B = o, q A bM a m b = e = ±1. (i) 

where Gab (ff/ii/) is the metric of the bulk (brane) space V m (V m -i) in arbitrary coordinates, {y A } 
({x^}) is the basis of the bulk (brane) and N denotes the vector normal to the brane. The perturba- 
tion of V m —i in a sufficiently small neighborhood of the brane along an arbitrary transverse direction 
£ is given by 

z A (x^,o = y A + (c ( y) A , (2) 

where C represents the Lie derivative. By choosing £ orthogonal to the brane, we ensure gauge 
independency [21] and have perturbations of the embedding along a single orthogonal extra direction 
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N giving local coordinates of the perturbed brane as 

^(^0=^ + £A^OO, (3) 

where £ is a small parameter along N A that parameterizes the extra noncompact dimensions. One 
can see from equation (2) that since the vectors N A depend only on the local coordinates x^, they 
do not propagate along the extra dimensions 

J\f A {x fl ) =M A + Z[N,N] A = N A . (4) 

The embedding equations of the perturbed hypersurface V m -i is given by 

g AB z A ll z B ^ = g flu , g AB z A ^ B = o, g AB N A N B = e = ±1. (5) 

Denoting by K^ v = ~gABy A ^ E the extrinsic curvature of V m -i, and using (3), the perturbed 
metric g^ u can be written as 

g»u = g^u - 2^ + fg a(3 K^K ul3 , (6) 

with the extrinsic curvature of the perturbed brane given by 

= -g AB Z A ^N B u = - iK m K%. (7) 

Comparison of equations (6) and (7) then leads to the York relation, describing the metric evolution 
with respect to the perturbation parameter £ 



K„ - - g-^-. (8) 



The components of the Riemann tensor of the bulk written in the the embedding {Z A , Af A } lead to 
the Gauss-Codazzi equations given by [24] 



TI AB cd2 a Z^Z^Z^ = R a /3j5 + ZK^K^s, (9) 



T^-abcd2 a M b Z^Zjg = Kj3y. a - K ar>f3 , (10) 

together with [19] 

Kabcd Z A tl Z c u M B M D = -C M K^ v + K^ a K a vl (11) 

where TZabcd and R a p^8 are the Riemann tensors for the bulk and the perturbed brane respectively 
and £/v denotes the lie derivative in the AA-direction. Using this projection, the m-dimensional 
Riemann curvature and its contractions (the Ricci tensor and scalar curvature) are described by the 
(m — l)-dimensional variables on the brane with the normal M A as [19] 

TIabcd = Rabcd - K A cK B d + K AD K B c - NaKbd,c + NaKbc-,d + N B K AD -c - N b Kca-,d 

-N C K BD ;A +NcK A D;B +N D K B C;A ~ M D K AC . B 

+N a M c KbeK e d - M a M d KbeK e c - N B N C K AE K E D + M B M D K AE K E C 
-NaNcCnKbd + NaNdCnKbc + N b N c CnK A d - MbM d C n K A c, (12) 

TZab = Rab - KKab + 2K AC K C B + Ma (k c b . c - K, B ) + N B (k c a . c - K, A ) 

+N A N B K CD K CD - C M K AB - M A M B g CD C M K CD , (13) 

U = R- K 2 + 3K CD K CD -2g CD C N K CD . (14) 

The Einstein-Gauss-Bonnet equation is quasi- linear [11], which means that apart from non-singular 
terms given by the (m — l)-dimensional variables, it contains only linear terms of Cj^K AB with no 
quadratic terms appearing. 
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3 Effective field equations 

We consider the total action for space-time (M,Gab) with boundary (S,^) as 

S =T~7 I d m X^G(TZ - 2AW + (3C GB ) + / d m - l x^(C snT{&cc + C m ). 
2a* Jm Js 

Variation of the total action gives our basic field equations as 

G% + K^G AB + pH% = a* Sab, 



where 



G 



n 



(b) 

AB 

(b) 
AB 



and 



TZ>ab - -^GabT^, 

TITIab — 2TZac'R- C b ~ 2TZ cd TZacbd + Ha DE 'R-bcde 



c G b = iz 2 - m AB n AB + KabcdTi abcd 



-^Gab^-gb, 



(15) 

(16) 

(17) 
(18) 

(19) 



is the Gausss-Bonnet term. In the above equation a* = M m--i {M* is the fundamental scale of energy 

in the bulk space), (3 is the Gauss-Bonnet coupling constant with dimension (length) 2 , is the 
cosmological constant of the bulk and Sab consists of two parts 



Sab = Tab + -^VGab, 



(20) 



where Tab = ~^j^£k + QAB^m is the energy-momentum tensor of the matter confined to the brane 
through the action of the confining potential V. We require V to satisfy three general conditions: 
firstly, it has a deep minimum on the non-perturbed brane, secondly, depends only on extra coordi- 
nates and thirdly, the gauge group representing the subgroup of the isometry group of the bulk space 
is preserved by it [16]. Substituting relations (12), (13) and (14) into Eq. (16), we find the effective 
equations on the brane as 

iV - \ p ^v + K W K% - g^K ap K^ - C M K^ U + g^g p(T C M K pa 



+ 



2(3 (H^ - PCmK^ + 2P/C N K pv + 2P v p L M K pp + W pu pa C N K pa ) 



a*S AB Z A ,Z B u ~ A {6) 



P + I3(P 2 - 4P a(3 P aP + P aPl5 P aM6 ) = -2a*S AB M A N B + 



(21) 
(22) 



where 



Pa(3fS 
PaP 
P 



H, 



pv 



W pa 

pV 



9 Pa Pa P Po = Ra/3 - KK af3 + K^K^, 

g a/3 P a p = R-K 2 + K a0 K 
PP 



2 r . T saf3^ 
2{PppP P v ~i~ P P PjjLpua) ~i~ PiipaK.Pi/ 

p ™ + 2K aP K a ^P tu/ + PK pp K" u 



- 2(K pp K? a Pl + K vp K p a P^) - 2K^K K °P ppv<7 



4^ 
9pu 



P 2 - 4P a(S P a P + P af s lS P a ^ 5 



Pgpvg pa - 2Pp U g pa - 2 9iXV P pa + 2P lwlvP g a Pg 



OLp j3a 



(23) 



(24) 
(25) 
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In order to find the effective equations on the brane, we have to replace the terms Cj^K pu in Eq. 
(21) with the (m — l)-dimensional variables on the brane. Using the decomposition of the Riemann 
tensor into the Weyl curvature, the Ricci tensor and the scalar curvature 

2 2 
Kabcd = Cabcd - ( m _ 2 ) (Sa[dKc\b - Gb[dK C ]a) - ( m _ x)( m _ ^ n (.^A[cGD]B), (26) 



we obtain 

C ^ K ^ ~ uT_ ^ 9^g aP C M K aP = - jg— + K^KP, - j-1 ^ \p pu - | 



p,v j -. \ PQ\iv 



{m-l)^ yv ap (m-3) ^ ' w v (m - 3) 

" (^TijVV". (27) 

where 

£„» = Cabcd N b N d Z a ^ Z c v , (28) 

is the electric part of the Weyl tensor Cabcd- However, since Eq. (27) is a trace free equation, we 
cannot fix Cj^K^ by it. We have to find g a @ Cj^K a p from the other independent equation. If we 
take the trace of our basic Eq. (16), we find 

(m - 2)11 + (3{m - 4)C GB = -2a*S + 2mA (b) . (29) 

Substituting Eqs. (9)-(ll) with Eq. (27) into Eq. (29), we obtain 

g C M K a ,--+K a ,K + — _ 2)+ — _ 4)] + 2 p _ 2) - pp{m _ m , (30) 

where 

I = P 2 - 8P a(S P a P + P a( s l5 P a ^ 5 - l2P pa £ pa . (31) 
From Eq. (27) with Eq. (30), we then obtain 

r K ( m ~ 2 ) c 1 ( p Ip A | (a*S-mA( fe )) 
^ ^ ~ (m - 3) ^ (m - 3) v ^ 2 ^7 + (m - 1) [(m - 2) + /?P(m - 4)] 5 ^ 
, K KP , (rn - 4)/?J 

+A W ,A I/ + _ x) _ 2) + ^ p(m _ 4)] 9ilv . 

Substituting Eq. (32) into Eq. (21), we obtain the effective gravitational equations on the brane as 

(m-2) ua(n™\n™\ 2(m-4) 2 /? 2 / / 1 \ 

( m _ 3 ) ^ + ^ + ^ r ^ + H ^ ) (m-l) [(m - 2) + /?P(m - 4)] ^ 4^7 

^ -a*S -AW, ,, 4(m-4) (a»g-mA(»>) / _1 \ 
a ^ + ^ ( m - 1) [( m - 2) + (3P(m - 4)] ^ 4^7 

9iiv , (33) 



2(m - 3) 
(a*S-mA^) 



(m-l) 
where 



#(i) _ 2 P »P°fo < m - 2 ) pp°p | 6 P p 4(m-2) P p 



_l_ g/^ 



( m _ 3 ) 1 (m-3) (m-3) 

(m 2 - 9m + 6)P 2 + 4(m 2 - 4m + 7)P a/3 P al3 



2(m- l)(m-3) 

- (m 2 - 5m + 6) P^ 7<5 P Q/375 ] , (34) 
(2) _ 4(m - 2) / > 



(m — 3) 

+ t 7T7 7p;9pvPpcr£ p H t 7prP£pv- (35) 
(m — l)(m — 3) (m — 3) 
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As was mentioned in the introduction, localization of matter on the brane is realized in this model 
by the action of a confining potential. Let us take 

aT ^ = ^a r ,„, A = -^A<»>, (36) 
(m — 2) (m — 1) 

where a is the scale of energy on the brane. Now, we rewrite Eq. (33) as an Einstein-type equation 
with "correction" terms. From Eqs. (22) and (33) we obtain 



[(m - 2) + (3P{m - 4)] 



(m 

Aim — 4) 4m(m — 4) , 

-y yCtT + — i -^A 

(m — 1) (m — 2) 



(37) 



where 



H-(i) - ^ ?)p fl P^7_ 4 pP^p , ? pp _Ap pP 

n n.v — (jn — 2) 1 ppua (m — 2) u w " 



9/j.v 



2(m-l)(m-2) 
+ (-2m 2 + 9m - 9)P a ^ 5 P a ^ S 



(-5m + 3)P 2 - 4(2m 2 - 8m + 10)P Q/3 P a/3 

2/? (m-4) 2 (m-3) 



[(m - 2) + /?P(m - 4)] (m - l)(m - 2) 
x [P^u ~ \P9^ (P 2 ~ 8P af3 P a P + P af3y5 P a ^ 5 ) , (38) 

H^J = —4 [p w £ p v + Pv P £ P n + Pfipi/ a £ p ) + _ ]_^ m _ r^9iivPpo£ p + 2P£ /iI , 



24/3 (m-4) 2 (m-3) / 1 , p 

[(m - 2) + /?P(m - 4)] (m - l)(m - 2) ~ 4^"" ' pa ' 



Here, A is the effective cosmological constant in four dimensions with Q^ u being a completely geo- 
metrical quantity given by 

Qtxv = KKp V - Kp p K p v + ^(K a pK a P - K 2 ) 9fiu . (40) 

A brief discussion on the energy-momentum conservation on the brane would be in order here. 
The contracted Bianchi identities in the bulk space G AB ^ = and Tt AB ^ = 0, using Eq. (16), 
imply 

+ 1 vg AB\ = Q (41) 
2 J -A 

Since the potential V has a minimum on the brane, the above conservation equation reduces to 

= 0, (42) 

As a check of our calculations so far, we note that Eq. (37) reduces to Eq. (32) given in [17] in a 
5-dimensional bulk with no Gauss-Bonnet term, that is (3 = 0. As can be seen from the definition 
of Qpu, it is independently a conserved quantity which, according to [20], may be interpreted as an 
energy-momentum tensor of a dark energy fluid representing the X-matter, also known as the 'X- 
Cold-Dark Matter' (XCDM). This matter has the most general form of the equation of state which 
is characterized by the following conditions [22]: first it violates the strong energy condition at the 
present epoch for uj x < —1/3 where p x = to x p x , second, it is locally stable i.e. c 2 = 5p x /Sp x > 0, and 
third, causality holds good, that is c s < 1. Ultimately, we have three different types of "matter" on 
the right hand side of Eq. (37), namely, ordinary confined conserved matter represented by t^ u , the 
matter represented by Q pu which will be discussed later and finally, the Weyl matter represented by 

£jXV- 
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4 Modified Friedmann equations 



In this section we will analyze the influence of the trace of t^, the extrinsic curvature terms and 
higher curvature terms on a FRW universe, regarded as a brane embedded in an m-dimensional bulk 
with a Gauss-Bonnet term. The FRW line element is given by 



ds 2 = -dt 2 + a{tf 



dr 2 



+ 



r 2 (d6 2 + sin 2 Ody 2 ) 



1 — kr 2 

The confined source is the perfect fluid given in co-moving coordinates by 
Tnu = {p + p)Uf,u u + pg^, Uf, = -<5°, p = (7 - l)p. 



(43) 



(44) 



Let us consider an AdS m , dS m or flat bulk, so that £^ u = 0. For late times this assumption seems 
reasonable because the effects of such a term is negligible. The Codazzi equations (10) reduce to 



K, 



k — o 



(45) 



Using the Yorks relation, it follows that in the FRW space-time (diagonal metric), K^ v is diagonal. 
After separating the spatial components, the Codazzi equations reduce to 



K, 



k r a 

v<J pp 



K 



p,p,v 



K T a 



(46) 



K^q - = -ad (<fj# + rXSl + r 2 sm 2 6ffi) K 00 . (47) 

The first equation gives Kn jV = for v 7^ 1, since K\\ does not depend on the spatial coordinates. 
After defining K\\ = b(t), where b(t) is an arbitrary function of t, the second equation gives 

Koo = - 1 A(-)- (48) 



a dt \a 

For /j,, v = 2, 3 we obtain K22 = b(t)r 2 and K33 = b(t)r 2 sin 2 9 and generally (//, v 7^ 0) 

K„v = Xgpv (49) 
cr 

We find from Eq. (40) that 



Q^^-^U^-bAg^, Qoo = ^- (50) 



Denoting h = | and H = ^, the components of become 



6 : 



» 2 



2^-l)^, Qoo = ^- (51) 



It would now be interesting to see how the above geometrical interpretation is compared with the 
X-matter explanation. To this end we define Q^ v as a perfect fluid and write 

Qpu = ^ \{p x + Px)upU u + p x g^ u ] , Px = {lx - l)Px- (52) 
Comparison with the components of given by Eq. (51) gives 



,.. = -14(4-,), * = < 53 > 
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Use of the above equations leads to an equation for b(t) 



» = 1(4 -,*(»))* (54) 

It is interesting to note that this equation resembles one of the phenomenological candidates for dark 
energy, the X-matter [22], but in our case this field has a fundamental geometrical justification for 
the equation of state, having been derived from the term Q„v in the Einstein equation (37), itself a 
result of the extrinsic curvature. If 7a; is taken as a constant, the solution for b(t) is 

b(t) = boa{t)^ 4 - 3 ^\ (55) 

where 6o is an integration constant. Using Eq. (53) and this solution, the energy density of XCDM 
becomes 

Px = ^a~^. (56) 
a 

It should be noted that we have considered an AdS m , dS m or flat bulk and consequently = 0, 
therefore Eq. (41) leads to H$ = 0. 

Now, using these relations and Eq. (37), the Friedmann equations, to first order in (3 become 

3H 2 + f{a) + (3 \H 2 A{a) + HB(a) + 18ii~ 4 + AH 2 + 1QHH 2 + f(a)} + 0{f3 2 ) = 0, (57) 



-2H - 3H 2 + g(a) + (3 
where 



H 2 C{a) + HD(a) - 18H 4 - ^H 2 - ^-HH 2 + g{a) 



+ O(/3 2 ) = 0, (58) 



A(a) 
B{a) 

f(a) 
Ha) 



and 



20k 



4A + — - 12&2 a -37,( 3 _ 2 lx ) - apoa _37 (37), 
1^ + !* _ 46ga- 3 ^(4 - 3 7 .) - 2apoa-^, 



Ah — » — 36na 3 ^ x — -apQa 



-3 7 



-(4 - 87x)A6ga-^« + 36ga"^ (+3 7 2 - 87. + 6) + 3ap ^ 3(7+7a;) (7 - 7*y 



+ 6 Iti 



3^ 



+ fcap (2 - 3 7 )a~ 37 " 2 - fc&§(20 - 12j x )a' 3 ^- 2 , 



52k 



C{a) = "4A-^2 +46 2 a- 3 ^(9-7 7 x)-apo7«" 37 , 

8A 40/e - /14 \ 2 ~ 
= - T -^ + 46 2 a" 3 -( y -2 7 ,)--ap a^ 

g(a) = -A-^+3b 2 a- 3 ^(l- 7x )-^poa- 3 \ 

g(a) = -bta'^ (e 7 2 - 28 lx + lg) + apo&ga" 3 ^) ( 7 - 7x ) + A^a"^- (4 - ^) 



(59) 



(60) 



Equations (57) and (58) now constitute our Fridmann equations, having been modified by a number 
of extra terms. Such terms may be used to offer an explanation for the X-matter discussed earlier. 
In the next section, we discuss the implications of such terms on the cosmology of our model. 
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5 Cosmological implications 



Recent observational data indicate that our universe is expanding with a positive acceleration. This 
acceleration is explained in terms of the so-called dark energy which could result from some exotic form 
of matter with negative pressure. The nature of such dark energy constitutes an open and tantalizing 
question connecting cosmology and particle physics. The simplest form of dark energy is the vacuum 
energy (cosmological constant). However, this scenario faces the two well known cosmological and 
coincidence problems [25]. Another possible form of dark energy is provided by scalar fields. Dark 
energy can be attributed to the dynamics of a scalar field, called the quintessence [26]. Also a 
phenomenological explanation based on current observational data is given by the X-matter model 
which is associated with an exotic fluid characterized by an equation of state p x = w x p x , where the 
parameter w x can be a constant or more generally a function of time [27]. In this section, following 
the method developed in [17], we show that within the context of the present model, it is possible to 
have a universe exhibiting accelerated expansion without the need to resort to a cosmological constant 
term or a scalar field or any other kind of dark energy. 

Let us proceed by taking equations (57) and (58) from which, the Hubble parameter can be 
written as 



H 2 



1 

210 

3 3 



2A _ ^ + 2Wb 2 a- 3 ^ - 3 (j - ^) ap^ 
1 + 2/3 



70/3 70/3 



-8A+ (^±l)ap a- 3 



(61) 



We note that for /3 — ► the negative sign yields a singular term in the above equation and, therefore, 
we select the positive sign. A qualitative classification of the solutions on the basis of different values 
of the parameter j x can be achieved without solving the equation. If we define the effective potential 



V{a) = k-bla-^ + ^(l- 2 ^)a-^ 



we then have 



+ ^-w{ 1 +i^ +7)Qwr3, } i - <62) 



a 1 + V(a) = 0. (63) 



Now, it is possible to deduce the qualitative behavior of the scale factor a(t), keeping in mind that 
d 2 should be positive. This behavior is much dependent on the range of the values that ^y x can take. 
We distinguish the following possibility for having an accelerating universe [23] 

< 7* < \- (64) 

We also note that to avoid imaginary values for the effective potential, (3 > 0. The behavior of the 
potential for k = 0, ±1, is illustrated in figure 1. This figure shows that, within the context of the 
present model, the universe exhibits accelerated expansion in the case of a vanishing cosmological 
constant with (3 > 0. This case is what we expect from the slope of the expansion in heterotic 
string theory. It should be mentioned that Gauss-Bonnet theories with a negative coupling constant 
have been suggested [28] in the past where cosmic acceleration can be achieved, but this is not well 
motivated by string theory and has problems with stability [29]. 

The above equation cannot be solved in closed form. However, in two extreme cases corresponding 
to small and large a(t), exact solutions can be easily found 

»W=(|) § (^)*rf. for small a, (65) 



9 




and 

b 2 

o(t) = -j-t 2 , for large a. (66) 

The first solution is of the Einstein de-Sitter type, while the second represents an evidently inflationary 
of the power-law type. This means that in our model the universe starts as decelerating and finally 
ends up as accelerating. In the simplest FRW cosmological models with a one-component fluid filling 
up the universe such behavior is not possible. In the next section we discuss the observational 
parameters of our model. 

6 Age of the universe 

The age of the universe in FRW models is given by 

t^if- rr, (67) 

HoJo {^ + (i-n m )Y 

where Hq 1 = 9.8 x 10 9 h _1 years and the dimensionless parameter h, according to modern data, is 
about 0.7. Hence, in the flat matter dominated universe with Qtotal = 1 the age of the universe would 
be only 9.3 Gyr, whereas the oldest globular clusters yield an age of about 12.5 with an uncertainly 
of 1.5 Gyr [30]. With the assumption of the cosmological constant, the WMAP data [31] yields 
to = 13.7 ± 0.2 Gyr, and the best estimate of the dynamical age of the universe, coming from CMB 
data, is t = U± 0.5 Gyr [32]. 

It would now be interesting to compare the above results and those predicted by the present work 
to the results obtained in [17] where the same mechanism for localization of matter is used in the 
absence of the Gauss-Bonnent term and the age of the universe is given by 

= w [ ; dx (68) 
HoJo 0+^*) 2 

for a flat, matter dominated universe with U m = 0.3 and U x = 0.7. This leads to a prediction for the 
age of the universe of about 12 Gyr. Here, with the GB term present, we introduce the dimensionless 
cosmological parameters 

Mm = oi j2 ' U * = "nTT2 > M GB = pH Q , (69) 

where the subscript zero refers to the present values of the cosmological quantities. The modified 
Friedmann equation (61) with A = and k = is then given by 



H 2 \oqJ 70 v 2 j \aoJ 70U GB S 



37 

1 + 2 (47 + 7) u GB n m ( — ' 



«0 



^ 1. 

2 



(70) 
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This equation leads to the constraint equation that should be satisfied amongst the cosmological 
parameters defined above when a = ciq, that is 



1 = Q x - — 
70 



(7 - ^pj n m - {1 - [1 + 2 (4 7 + 7) n GB n m }^} , (71) 



which reduces the number of independent Q parameters. We therefore find the age of the universe 
by direct integration of the Friedmann Eq. (70) 

= 7T f " —™ 

" J0 {sixx- 3 T+ 2 - 1,(7 - S2)a m x- : <-'+' - „gj [1 - (i + 2(4 7 + 7)n GB n m x-^] } 2 

The results are summarized in table 1. It shows that the age of the universe in our model is longer 
than the FRW model and brane models without the Gauss-Bonnet term in the bulk action [17]. 
Noting that the dynamical age of the universe depends on the rate of the expansion, we realize that 
in the presence of the GB term, the universe accelerates faster than brane models without the GB 
term. We have plotted the age of the universe in three models as a function of the energy density 
parameter Q. m in figure 2. Figure 3 shows the age of the universe as a function of Q. m , corresponding 
to the values of ^ x and 7. In figure 4 , we have plotted the age as a function of two parameters, 
( r y x ^m) and (7,O m ). Note that the age of the universe increases with decreasing 7x , 7 and il m . 



72 t Q H Q t (Gyr) 



0.1 


1.061 


14.54 


0.2 


1.037 


14.21 


0.3 


1.010 


13.83 


0.4 


0.978 


13.41 


0.5 


0.942 


12.91 


0.6 


0.901 


12.34 



Table 1: Age of the universe for various values of 7; 



At this point, it would be appropriate to compare these results with other brane world models 
where the Gauss-Bonnet term is absent. A very noticeable and popular scenario in this connection is 
the Randall-Sundrum brane model where the effects of the brane parameters and dark energy on the 
age of the universe have been studied. It has been shown that the effect of the quadratic term p 2 , 
resulting from the imposition of the Israel junction conditions, in the energy density term is to lower 
significantly the age of the universe. The quadratic term contribution in the energy density even for a 
small negative pressure contributes effectively as the positive pressure, and makes brane models less 
accelerating. This problem can be avoided if we accept the dark energy p = — |p (phantom matter) 
on the brane, since it has a very strong influence on increasing the age [33]. 

Another scenario which has attracted a considerable amount of attention in recent years is the 
so-called DGP model, proposed by Dvali, Gabadadze and Porrati (DGP) [34] and generalized to 
cosmology by Defayet [35] . This proposal explains the observed late time acceleration of the expansion 
of the universe through a large scale modification of gravity coming from the "leakage" of gravity 
at large scales into an extra dimension, without requiring a non-vanishing cosmological constant 
[36]. In other words, the bulk gravity sees its own curvature term on the brane as the cosmological 
constant and accelerates the universe. From the observational viewpoint, it has been shown that such 
models are in agreement with the most recent cosmological observations. For example, in [36], the 
authors show that constraints from the SNe Ia+CMB data require a flat universe with Q m = 0.3 and 
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0.2 0.4 
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Figure 2: The age of universe with the GB term (dot-dashed line), without the GB term (dashed line) and as predicted 
by the FRW model (solid line), as a function of Q m torj = 1, 7 X = ~ and Q,gb = 10 -3 . 





0.2 0.4 



Figure 3: Left, the age of universe as predicted by the present model as a function of fl m for different values of 
j x = 0.5, 0.3, 0.1, bottom, middle and top curves respectively with 7 = 1 and Q.gb = 10 -3 and right, the age of universe 
as a function of Q. m for different values of 7 = 2, 1.5, 1, bottom, middle and top curves respectively with 7^ = | and 
SIgb = KT 3 . 

£l rc = 0.12, where Vt r< , is the density parameter associated with the crossover distance between the AD 
an 5-D gravities. Also, the age of the universe in the context of these models have been studied and 
shown that for a fixed value of fi m , the predicated age is longer for larger values of f2 rc and therefore 
these models are more efficient in addressing the problem of estimating the age of the universe [37] . 

7 Conclusions 

In this paper, we have studied the Gauss-Bonnet brane world cosmology where the matter is confined 
to the brane through the action of a confining potential, rendering the use of any junction condition 
redundant. The modified Friedmann equations on the brane were obtained and shown to be modified 
by various terms. As a notable consequence, it was shown that dark energy can emerge as a result 
of the extrinsic curvature in the Gauss-Bonnet brane world. Therefore, this scenario features the 





Figure 4: Age of the universe as predicted by the present model as a function of -y x and fl m , left and as a function of 
7 and f2 m , right. We have taken Qgb = 10~ 3 . 
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possibility of an accelerated expansion at the late stage of the cosmic evolution without the need to 
invoke either a cosmological constant or a quintessence component. It turns out that if the Gauss- 
Bonnet coupling is positive, the universe undergoes an accelerated expansion. The addition of the 
Gauss-Bonnet term has resulted in an estimate for the age of the universe that is compatible with 
modern observational data, improving on our estimates in a previous work where the same mechanism 
was used to localized the matter on the brane in the absence of the Gauss-Bonnet term. 



References 

[1] L. Randall and R. Sundrum, Phys. Rev. Lett. 83 (1999) 3370 [hep-ph/9905221]. 

[2] L. Randall and R. Sundrum, Phys. Rev. Lett. 83 (1999) 4690 [hep-th/9906064] . 

[3] P. Brax and C. van de Bruck, Class. Quant. Grav. 20 (2003) R201 [hep-th/0303095] , 

R. Maartens, Reference Frames and Gravitomagnetism, ed. J Pascual- Sanchez et. al. (World 
Scientific,2001), p.93-119 [gr-qc/0101059], 

D. Langlois, Prog. Theor. Phys. Suppl. 148 (2003) 181 [hep-th/0209261]. 

[4] B. Zwiebach, Phys. Lett. B 156 (1985) 315, 

D. G. Boulware and S. Deser, Phys. Rev. Lett. 55 (1985) 2656 . 

[5] A. H. Chamseddine, Phys. Lett. B 233 (1989) 291, 

F. Muller-Hoissen, Nucl. Phys. B 346 (1990) 235. 

[6] N. D. Birrell and P. C. W. Davies, Quantum Fields in Curved Space (Cambridge University 
Press, Cambridge, England, 1982). 

[7] D. Lovelock, J. Math. Phys. 12 (1971) 498, 

M. Farhoudi, Gen. Pel. Grav. 38 (2006) 1261 [physics/0509210]. 

[8] J. E. Kim and H. M. Lee, Nucl. Phys. B 602 (2001) 346 , ibidB 619 (2001) 763 [hep-th/0010093]. 

[9] K. A. Meissner and M. Olechowski, Phys. Rev. Lett. 86 (2001) 3708 [hep-th/0009122], ibid Phys. 
Rev. D 65 (2002) 064017 [hep-th/0106203]. 

[10] Y. M. Cho and I. P. Neupane, Int. J. Mod. Phys. A 18 (2003) 2703 [hep-th/0112227] , 

Y. M. Cho, I. P. Neupane and P. S. Wesson, Nucl. Phys. B 621 (2002) 388 [hep-th/0104227], 

N. E. Mavromatos and J. Rizos, Phys. Rev. D 62 (2000) 124004 [hep-th/0008074], 

I. P. Neupane, JEEP 0009 (2001) 040 [hep-th/0008190], ibid Phys. Lett. B 512 (2001) 137 

[hep-th/0104226]. 

[11] N. Deruelle and M. Sasaki, Prog. Theor. Phys. 110 (2003) 441 [gr-qc/0306032]. 

[12] B. Abdesselam and N. Mohammedi, Phys. Rev. D 65 (2002) 084018 [hep-th/0110143], 
N. Deruelle and T. Dolezel, Phys. Rev. D 62 (2000) 103502 [gr-qc/0004021], 
S. Nojiri and S. D. Odintsov, JEEP 0007 (2000) 049 [hep-th/0006232], 
J. E. Lidsey, S. Nojiri and S. D. Odintsov, JEEP 0206 (2002) 026 [hep-th/0202198], 
C. Charmousis and J. Dufaux, Class. Quant. Grav. 19 (2002) 4671 [hep-th/0202107] , 

C. Germani and C. Sopuerta, Phys. Rev. Lett. 88 (2002) 231101 [hep-th/0202060], 

G. Kofinas, R. Maartens, and E. Papantonopoulos, JEEP 0310 (2003) 066 [hep-th/0307138], 

D. Iakubovskyi and Y. Shtanov, Class. Quant. Grav. 22 (2005) 2415 [gr-qc/0408093]. 

[13] W. Israel, Nuovo Cim. B 44 (1966) 1, 

S. C. Davis, Phys. Rev. D 67 (2003) 024030 [hep-th/0208205], 

E. Gravanis and S. Willison, Phys. Lett. B 562 (2003) 118 [hep-th/0209076]. 



13 



[14] R. A. Battye and B. Carter, Phys. Lett. B 509 (2001) 331 [hep-th/0101061]. 

[15] V. A. Rubakov and M. E. Shaposhnikov, Phys. Lett. B 125, (1983) 136. 

[16] S. Jalazadeh and H. R. Sepangi, Class. Quant. Grav. 22 (2005) 2035 [gr-qc/0408004] . 

[17] M. Heydari-fard, M. Shirazi, S. Jalazadeh and H. R. Sepangi, Phys. Lett. B 640, (2006) 1 [gr- 
qc/0607067]. 

[18] T. Shiromizu, K. Maeda and M. Sasaki, Phys. Rev. D 62 (2000) 024012 [gr-qc/9910076]. 
[19] K. Maeda, T. Torii, Phys. Rev. D 69 (2004) 024002 [hep-th/0309152]. 

[20] M. D. Maia, E. M. Monte, J. M. F. Maia and J. S. Alcaniz, Class. Quant. Grav. 22 (2005) 1623 
[astro-ph/0403072]. 

[21] M. D. Maia and E. M. Monte, Phys. Lett. A 297 (2002) 9 [hep-th/0110088]. 

[22] M. D. Maia, E. M. Monte, J. M. F. Maia, Phys. Lett. B 585 (2004) 11 [astro-ph/0208223]. 

[23] M. Turner and M. White, Phys. Rev. D 56 (1997) 4439 [astro-ph/9701138], 

T. Chiba, N. Sugiyama and T. Nakamura, Mon. Not. R. Astron. Soc. 289 (1997) L5 [astro- 
ph/9704199]. 

[24] L. P. Eisenhart 1966 Riemannian Geometry, (Princeton University Press). 

[25] P. J. E. Peebles and B. Ratra, Rev. Mod. Phys. 75 (2003) 559 [astro-ph/0207347], 
T. Padmanabhan, Phys. Rept. 380 (2003) 235 [hep-th/02 12290]. 

[26] R. R. Caldwell, R. Dave and P. J. Steinhardt, Phys. Rev. Lett. 80 (1998) 1582 [astro-ph/9708069], 
A. R. Liddle and R. J. Scherrer, Phys. Rev. D 59 (1998) 023509 [astro-ph/9809272]. 

[27] Z. H. Zhu, M. K. Fujimoto and D. Tatsumi, Astron. Atrophys. 372 (2001) 377 [astro-ph/0107234], 
P. S. Corasaniti, M. Kunz, D. Parkinson, E. J. Copeland and B. A. Bassett, Phys. Rev. D 70 
(2004) 083006 [astro-ph/0406608]. 

[28] M. H. Dehghani, Phys. Rev. D 70 (2004) 064009 [hep-th/0404118]. 

[29] D. G. Boulware and S. Deser, Phys. Rev. Lett. 55 (1985) 2656. 

[30] L. M. Krauss and B. Chaboyer, Science 299 (2003) 65 [astro-ph/01 11597]. 

[31] W. L. Freedman and M.S. Turner, Rev. Mod. Phys 75 (2003) 1433 [astro-ph/0308418]. 

[32] L. Knox, N. Christensen and C. Skordis, Astrophys. J. 563 (2001) L95 [astro-ph/0109232]. 

[33] W. Godlowski and M. Szydlowski, Gen. Rel. Grav. 36 (2004) 767 [astro-ph/0404299], 

M. P. Dabrowski, W. Godlowski and M. Szydlowski, Mod. Phys. D 13 (2004) 1669 [astro- 
ph/0210156]. 

[34] G. Dvali, G. Gabadadze and M. Porrati, Phys. Lett. B 485 (2000) 208 [hep-th/0005016], 
G. Dvali and G. Gabadadze, Phys. Rev. D 63 (2001) 065007 [hep-th/0008054]. 

[35] C. Deffayet, Phys. Lett. B 502 (2001) 199 [hep-th/0010186]. 

[36] C. Deffayet, G. Dvali and G. Gabadadze Phys. Rev. D 65 (2002) 044023 [astro-ph/0105068]. 

[37] J. S. Alcaniz, D. Jain and A. Dev, Phys. Rev. D 66 (2002) 067301 [astro-ph/0206448], 
J. S. Alcaniz, Phys. Rev. D 65 (2002) 123514 [astro-ph/0202492]. 



14 



